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1 Introduction 

Inflation has been very successful in explaining several puzzles of the standard big bang 
scenario. But the most important success of the inflationary theory, is its prediction of 
almost scale invariant power spectrum of primordial fluctuations [1-5]. In the inflationary 
scenario the temperature fluctuations of cosmic microwave background (CMB) and the large 
scale structure (LSS) are directly related to the curvature perturbations produced during 
inflation. Current observations of CMB strongly support the presence of an almost scale- 
invariant power spectrum. 

From a theoretical point of view, inflation is important because it gives an opportunity 
to test predictions of quantum field theory in a curved space time. In this framework, the 
choice of vacuum for quantum fluctuations is ambiguous. Hamiltonian for the primordial 
fluctuations is time-dependent and hence the meaning of vacuum is not very clear. Power 
spectrum of primordial fluctuations is generally computed assuming that the fluctuations 
are initially in the Bunch-Davies vacuum state[6, 7]. Bunch-Davies vacuum state is the 
minimum energy eigenstate of the Hamiltonian in the infinite past and it is a reasonable 
choice as an initial state but not unique. It is somewhat of a philosophical question whether 
initial conditions are integral part of a theory or should be analyzed separately. There has 
been a great deal of work focused on modifications of Bunch-Davies vacuum state[8-19]. In 
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this paper, for simplicity we will assume that the Bunch-Davies state represents an empty 
vacuum state in the beginning of inflation. Even with that, there is no good reason why the 
initial state has to be a vacuum state. We do not know anything about the physics before 
inflation and a priori any excited state is as good an initial state as the vacuum state. In 
fact excited initial states can open interesting possibilities. We also believe understanding 
the initial state of the primordial fluctuations will be the first step towards the physics before 
inflation. 

In spite of the great success of the inflationary theory, it is always important to verify the 
validity of different assumptions. The purpose of this paper is to discuss the effect of relaxing 
the assumption of Bunch-Davies vacuum by choosing a general initial state built over the 
Bunch-Davies vacuum state. The effect of having thermal initial state on the power spectrum 
has already been considered in [20, 21]} In [23] authors have identified a family of excited 
states that are indistinguishable from the Bunch-Davies vacuum at the level of two-point 
function or three-point function. Similar work has also been carried out by Agullo and Parker 
in [24] but our approach in this paper is different. Today our universe is homogeneous and 
isotropic on the large scale and the power spectrum of primordial fluctuations is nearly scale- 
invariant from the size of the observable universe down to the scales of around a Mpc; as the 
authors show in [25]. Also the measurements of the bispectrum (three-point function) from 
the CMB[26] and LSS[27] indicate that the primordial fluctuations are nearly gaussian. The 
current observational constraint on the non-gaussianity parameter /tvl is very weak and 
remains the best constrained non-g aussianity parameter: -5 < < 59 (WMAP7+SDSS). 
Quasi-de Sitter inflation (or slow roll inflation) with Bunch-Davies vacuum state as the initial 
state for the fluctuations produces an almost scale-invariant power spectrum (n^ ~ 1) and 
a negligible non-gaussianity (/^£ ~ 0.02) [28]. As a first step, we will compute the power 
spectrum and the bispectrum with coherent state as the initial state and show that a large 
class of coherent states can produce exactly the same power spectrum and bispectrum as the 
vacuum state. 

We will then compute the power spectrum for quasi de Sitter inflation with a gen- 
eral initial state and find the constraint on the initial state by demanding scale invariance. 
Constraints on the initial state from current measurements of power spectrum are relatively 
weak. For quasi-de Sitter inflation, Bunch-Davies vacuum and coherent states are not the 
only states, in fact a large number of states are consistent with the observations. However, 
we will demonstrate that renormalizability of the energy-momentum tensor of the primor- 
dial fluctuations imposes some restrictions on the initial state. We will also argue that 
fNL ^ Ti^^ ~ '^'') ^olds even for a general initial state. Therefore, if a state produces an 
almost scale- invariant power spectrum then it can not produce large fj^l- 

The rest of the paper is organized as follows. In section 2, we review the quantization 
of the fluctuations in inflationary universe. In section 3, we compute the power-spectrum 
and the bispectrum with coherent state. Section 4 is devoted to the introduction of excited 
initial states and constraints on the initial state from observations. In section 5, we discuss 
renormalizability of the energy-momentum tensor of the fluctuations and constraints that 
imposes on the initial state. We end with a discussion on bispectrum in squeezed limit with 
a general state in section 6 and concluding remarks in section 7. 



^See [22] for a discussion on graviton spectrum from a thermal initial state. 
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2 Review of quantization of the fluctuations in inflationary universe 
2.1 Scalar field in FRW universe 

We start with the Lagrangian of gravity and minimally coupled real scalar field (for a nice 
review see [29]) 



S 



1 



-gd'^x 



1 



SvrG 



R - g'"'df,4>d,^ - 2V{(t>) 



A homogeneous background solution has the form 

ds'^ = -dt^ + a^{t)dx^ 
with a background scalar field (j){x,t) = (j){t). This background obeys the equations 



3H' 



+ 3H^ + V'{ 



02 + V{^) 



0, 



(2.1) 
(2.2) 

(2.3) 
(2.4) 



where H is the Hubble parameter H = d/a. Next we consider perturbations around the 
homogeneous background solutions. 



and the perturbed metric with only scalar perturbations is 



ds' 



[1 + 2<^)dt'^ + 2a{t){diB)dx'dt + a^it)[{l - 2^)6ij + 2dij E]dx' dx^ . 



(2.5) 



(2.6) 



We can avoid fictitious gauge modes by introducing gauge-invariant variables[30, 31]. One 
such variable is the comoving curvature perturbation 



71 = -^ + 



H . 



(2.7) 



Expanding the action (2.1), we get the second order action for the scalar fluctuations in terms 
of the gauge-invariant variable IZ 



So 



d^xa^ 



TZ^-a-^idiTZf 



(2. 



The next step is to quantize this system. Before doing that to simplify the action, we 
introduce Mukhanov variable v = zTZ and conformal time r, where z = a(j)/H. In terms of 
the equation of state parameter, z"^ = ^qOl^{^ + w)? That leads to the action 



S2 = \l drd^x [{v'f - {d,vf + vH{t)] , 



(2.9) 



where (...)' = dr{---) and 6(t) = [z" / z). We can define the Fourier transform of the v field in 
the standard way, 

t;(x,r)= / — — ^k(r)e^'^-^ (2.10) 



J (2vr)3 

^Equation of state parameter is defined hy p — wp, where p is the pressure and p is the energy density. 
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with v^{t) = T;_k(T) (because field v is real). And in terms of the Fourier transform the 
action (2.8) now reads 



S^ = lj^l drd'k [<(r)<(T) - k\Ur)vlir) + 5(r)^;k(r)<(r) 

From this action we get the following equation [5, 32] for 

v'^{t)+uI{t)v^{t) = 



(2.11) 



(2.12) 



with <^^(t) = k'^ — 6(t). If Uk{T) and uI.{t) are linearly independent and form a ba- 
sis in the space of complex solutions of equation (2.12) then the Wronskian P^[nA;,ti|,] = 
2iIm[u'f,{T)ul.{T)] / 0. Complex solution life(r) is called the mode function. Wronskian W is 
time- independent so we can always normalize the mode function Uk{T) by the condition 



/mK,(T)<(r)] = 1. 
The general solution of equation (2.12) can be written as 



(2.13) 



(2.14) 



where a^^ and a~^^^ are independent of r and aj^ = (oj^)* (because t;(x, r) is real). 
2.2 Quantization of fluctuations 

Field f (x, r) can be quantized in the standard way, just like harmonic oscillator (we will work 
in the Heisenberg picture). We introduce the commutation relation 



[{)(x, r), 7r(y, r)] = i5^(x - y), 

where tt = "0' is the canonical momentum. Now the equation (2.14) becomes 

1 



V2 



(2.15) 



(2.16) 



Now comoving curvature perturbation is an operator TZ = v/z. The commutation relation 
(2.15) gives commutation relations between dj^ and dk 



«ki , oL 



(27r)353(ki-k2), 



[dkijdka] = 0. 



The Hamiltonian of the system is 
1 1 



H{r) 



4 (27r)3 



dkd_k-Ffc(T) + d|^d^,^Ffc(r) + (^dkdj^ + dj^dkj -E'fc(r) 
Fkir) = {u'kf + ujIuI, Ek{T) = + ul\ul\. 



where, 



(2.17) 



(2.18) 



(2.19) 
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2.3 Bunch-Davies vacuum 

The next step is to define a vacuum state and find out the mode-functions that describe the 
state. The Hamiltonian explicitly depends on the conformal time r. So it is not possible 
to define a vacuum in a time-independent way. First we define the vacuum by the standard 
condition: for all k 

ak|0) = 0. (2.20) 

But this is not enough to specify the mode-functions. It is not possible to find a time- 
independent eigenstate of the Hamiltonian. So we take a particular moment r = tq, and 
define the vacuum as the lowest-energy eigenstate of the instantaneous Hamiltonian of the 
fluctuations at r = tq (we can always do that as long as w^(ro) > 0). That gives us following 
initial conditions for the mode functions 

4(ro) = ±iV^JJ^)e''^''\ n,(To) = ± , \ e^^W, (2.21) 

where X{k) is some arbitrary time independent function of k. In the limit when tq represents 
infinite past (i.e. tq — )• — oo), the vacuum is called Bunch-Davies vacuum state. In this limit, 
ujI = k'^ > and we can define vacuum by equation(2.20) for all modes. 

2.4 Power-spectrum for inflation 

During inflation,'^ we have w ^ —1 (see [33] for constraints on w) and 6(t) = (2/r^). Solving 
equation (2.12) with normalization condition (2.13) and initial conditions (2.21) (with the + 
sign and tq — ?■ — oo limit), we get 

ikr / ^ \ 

With this mode functions we can now compute the power-spectrum. For that we need to 
calculate following quantity 

(*k(T)l)k'(r)) = (0|t)k(r)7)k'(r)|0) = l{2TTf6'{k + k')\uk{T)\\ (2.25) 

where we got the last equation using (2.16). Before we compute (7^k(''")'7^k' (''")); should 
introduce some standard quantities 

(^k(T)7^kKr)) = (27^)3J3(k + k')P7^, A| = ^Pt^, n, - 1 = (2.26) 

where is called the scalar spectral index or tilt. Using equations (2.24-2.26), in the 
superhorizon limit (|A:r| <C l),we have 

n, = 1, Al = (2.27) 



3 



Here we will assume that V{(ji) is approximately constant and slow roll parameters |e|, 1, where 

^-^^ = ^(- + ^)-16^(y)'' (2.22) 
1 V" 



It is easy to check that H and are also approximately constants. 
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Where H is the Hubble parameter during inflation. Observation of CMB and LSS tehs us 
conclusively that the power spectrum of the fluctuations produced during inflation is almost 
scale-invariant (i.e. ~ 1). Here we have ignored the slow roll parameters because deviation 
from perfect scale invariance is small. In the flrst order in slow roll parameters, spectral index 
is given by 

ns = l- 6e + 2r]. (2.28) 



3 Coherent states 



In this section, we will compute the power spectrum and the bispectrum with a non- vacuum 
initial state namely, coherent state. Coherent state is a special state because it closely 
resembles classical harmonic oscillation. But the obvious question is: why coherent state? 
Because we do not know anything about the physics before inflation, a priori any excited 
state is as good an initial state as the vacuum state. In particular, coherent state is an 
interesting example of a non-vacuum state and we will show that a large number of coherent 
states are indistinguishable from the Bunch-Davies vacuum state. And it is a good first step 
towards a better understanding of the initial condition of primordial fluctuations. 

We can use aj^ operators to build excited states over the Bunch-Davies vacuum state 

The coherent state is deflned in the usual way, 

ak|C) = C{k)\C) (3.2) 
with (C|C) = 1. Now the vacuum state is the special case: C(k) = 0. 
3.1 Power-spectrum with coherent states 

Before marching onwards with our calculations, we want to impose the following constraint 
on the initial state. We want {TZ\i{T)) = 0, in superhorizon limit. That leads to the condition 

C*(-k) = C7(k). (3.3) 

It is important to note that renormalizability of energy-momentum tensor of the fluctuations 
(see section 5) constraints C(k) further: C(k) goes to zero faster than for large k. 
Now we can compute the power spectrum 

(^k(r)^k'(^)) = (C|^k(r)7tk'(r)|C), (3.4) 

using 

^k(r) = [aMr) + al^JZkir)] (3.5) 

where, TZkir) = (^^^ (l + With the help of equations(3.2,3.3) and commutation 
relations (2.17) in the superhorizon limit (|A;t|, |A;'r| <^ 1), we get the following result 

(7tk(r)7tkKT)) = h27rfJ^6\k + k'). (3.6) 
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Therefore, spectral index = 1.^ Note that equation(3.6) is exactly the same as the power 
spectrum with Bunch-Davies vacuum state. It is not very difficult to understand why we 
have the same result. One can think of coherent state as the zero-point quantum fluctuations 
around some classical state. So, the field '^^''^(t) in the coherent state can be written as 
7?.^°'*(r) = TI'^\t) + 'R,^'^{t), where classical part TZ'^\t) is obviously the expectation value 
{'R,^^{t)). And 'R,^'^{t) is the original quantum field but now in the vacuum state. Because 
of (3.3), it is very easy to check that the classical part TZ'^\t) decays and in the superhorizon 
limit W\t) vanishes. So, finally we end up with (7^^f (r)7t^°^(r)) = (Ttjif (t)'^™^(t)). 

3.2 Primordial non-gaussianities from coherent state 

Non-gaussianity (for a good review see [34]) is emerging as a very important tool for discrim- 
inating between different models of inflation. Even for slow roll inflation three-point function 
may contain more information about the initial state. As we have seen already that coherent 
state produces scale- invariant power spectrum just like the vacuum state. In this section we 
will compute the three-point function with coherent state (with slow-roll approximations) to 
see if there is any observable signature. To do that we will follow an approach similar to [24] 
by Agullo and Parker. At leading order in the slow-roll parameters, the third order action 
for the scalar fluctuations is given by [28] 

53 = SttG j d?xdT a\T) i -| j HTZ'^d-^TZ'^, (3.7) 



where, TZc is the redefined field 



1 _ . _o 1 



7^ = 7^e + ^ (3e - 2rj) TZl + -e 8^^ {11,8^11^) . (3.8) 

In momentum space the last equation becomes 

1 , f d^p 1 /■ d^p (k - p)2 , ^ 

= ^c,k + 4 (3e - 2?]) J ^7ec,p7ec,k-p + 2^ y (2vr)3 ^c.p7^c,k-p• (3.9) 

The interaction Hamiltonian can be found from S3 = — J drHint- In momentum space Hint 
is given by 

HUr) = ^«'(^) [Pj H I d'p^d'p2d'ps n'^^iT)n'^^{T)TZ'^.^iT)6\p,+p2+P3). 

(3.10) 

We will use time dependent perturbation theory to compute the momentum space three point 
function. But before that it is important to note that in the presence of interaction (3.10), 
{C\'7Z]i{T — )• 0)|C) 7^ 0. So the quantity we are interested in is not (C|7^ki(''")'^k2(''")'^k3(''")|C') 
but (C|7t^';^(r)7t£^(r)7t£^(T)|C), where 7tf^(r') = n^T')- {C\n^{T ^ 0)\C). Therefore, 
at leading order in slow-roll parameters,^ 

(C|<^(T)<^(r)7tg^(r)|C) =(C|7tk,(r)7tk,(T)7tk3(r)|C) 

Vl^kl(T)|C)(C|7^k2(r)7^k3(r)|C) + cyclic perm) . 

(3.11) 



*If we don't neglect the slow roll parameters then the spectral index, just like equation (2.28), is given by 
= 1 — 6e + 277. 

^ (C|7^k(r)|C) = for the Bunch-Davies vacuum and hence 7?.^''''(r') = 7tk(r')- 
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Using time-dependent perturbation theory we can calculate the three-point function (see 
Appendix A for explicit calculations) and in the limit r — t- 0, the final answer is 



(C|7l£^(r)7t£^(r)7t£^(T)|C) = {Infd^ik, + + ks) PR{k2) Pnih) 



2\^' + ' kl )^{kl+k2 + k3) kl 



+ cyclic perm. 
(3.12) 



Surprisingly, the last result does not contain the function C(k) and the result is exactly the 
same as the result for the vacuum state 

(C|7t£^(T)7t£f (r)^£^(r)|C) = (0|7tk,(r)7tk,(r)7tk3(r)|0). (3.13) 

So, it is not possible to differentiate between the vacuum state and the coherent states (with 
constraint (3.3)) from the power spectrum or the bispectrum. 

To understand the result (3.13), we can try to rederive it using a simple argument. Just like 
the two-point function case, we can decompose the field '^^y''^°^[j) in the coherent state 
into 7tf^''=°'*(r) = 7^f^''='(r) +^r(^)' ^ith 7^f^'''=V) = (Ttf ^''°^r)). By construction 
(7t^'^^'™^(r)) vanishes in the limit r — t- and as a result, in this limit TZ^^''^\t) does not con- 
tribute. Therefore, we have (^P^?/;c°'^(^)^P^2/;'=°'^(^)7^gy;c°^(^)) = ^7^mc(^)^j;ac(^^7^j;ac(^)^^ 
However, unlike the two-point function case, here it is not at all obvious that in the presence 
of interactions the field iz^y''^"^{r) can be written as '^^'^^•™'*(r) = W^'^ij) + TZ^^''^\t). 

4 General state 

Now as the next step, we will compute the power-spectrum with a general initial state and 
try to find out constraints on the initial state from observations. It is important to note that 
we are in the Heisenberg picture where states are time-independent. We can write down the 
most general state using equation (3.1) 

\G) = Y,C^\i^)- (4.1) 

v> 

And with this initial state we can compute the power spectrum lJZ\^{T)iZ\f^i {t)) = -£i {v]^{T)v\f^i {t)) 
using 

(vk(r)?Jk'(r)) = ^^^1^^^^ . (4.2) 

Therefore, 

(^k(r)i)k'(r)) = ^(27r)353(k + k')|txfc(r)|2 + yl(k,k')«' +^*(-k, -k!)ukUk' 

+B{-k, ]^)ukul, + B{-k', k)uluk' , (4.3) 



where. 
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Introducing k^, = \^kk', k = k + k' and Ak = k — k' and using equation (2.24), we can write 
(ik(r)*(r)) = \(2^mk + k')l (l + ^) + A(k,k')e-'«'i (l - - ^) 

In the superhorizon limit (|A;t|, |A;'r| ^ 1), we have, 

+ [-A(k, k') - A*(-k, -k') + B{-k, k') + i?(-k', k)] i^j^ + -^j^) +■■■ 

(4.6) 

where the dots indicate terms of higher order. Therefore, in the superhorizon hmit, at the 
leading order we have 

(7tk(r)7tkKT)) « l{27rfJ^6^ik + k') + J^[- ^(k, k') - ^*(-k, -k') 

2 (/)2/fc3 02^3 

+ B{-k,k') + B{-k',k)]. (4.7) 

4.1 Constraints from observations 

In this section, we will try to find the constraints on the initial state from observations. Our 
universe as we see it today, is homogeneous and isotropic on large scale. Temperature fluctu- 
ations of CMB and LSS are directly related to the curvature perturbations produced during 
inflation. Also the present observations of the CMB temperature inhomogeneities indicates 
the presence of almost scale-invariant spectrum of curvature perturbations. In general, the 
one-point function (7^k('?")) 7^ 0, even though TZk{T) is a perturbation.^ Homogeneity de- 
mands that (T^k (■'")) = but observation requires homogeneity only for the superhorizon 
modes. Therefore we will impose {T^kir)) = only in the superhorizon limit. From that we 
can write down the flrst constraint 

{G\al^\G) = (G|ak|G). (4.8) 

Homogeneity also demands that, in the superhorizon limit 

(7tk(r)7tk'(T)) = P{k, T)(53(k + k'). (4.9) 

Where P{k, r) is some arbitrary function of k and r. And the isotropy condition tells us 
that P(k, r) = P{k,T). Now demanding scale-invariance we get 

Pik,r) = {2nf^^, (4.10) 

®It is easy to understand this by considering the simpler case S — J dt [^x^ — V{x)'j. Action for the 

perturbation around the classical solution x is given by, S — J dt ^-^Sx — ^V" {x){Sx)^'^ . This perturbations 

can be quantized just like harmonic oscillator and thus in general {'^\Sx{t)\'ii) ^ for a general quantum state 
I*). 
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where W is some dimensionless constant independent of k. Factors (27r)'^ and H^/cfP' are 
included for later convenience. Comparing the last equation with the leading order term 
of equation (4.7), we also find that W does not depend on r and hence {T^k{T)TZ]i' (t)) is 
time- independent. Now comparing equation (4.9) with the equation (4.7), we get the second 
constraint equation 

- A{k, k') - A*(-k, -k') + 5(-k, k') + B{-k', k) = {2TrfW'5^{k + k'), (4.11) 
where W = {W — ^) is also a constant. And the power spectrum is given by 

(7^k(r)7^kKr)) = (2^)'^ + ^') '^'(k + kO, (4-12) 

where W is defined by equation (4.11). Equations (4.8, 4.11) are two constraints on the 
initial state. ^ It is reasonable to assume that the initial state |G) does not contain any 
excitation of infinitely large momentum. That is W = 0. In that case constraint (4.11) 
reads 

-^(k,k') -^*(-k,-k') + -B(-k,k') + -B(-k',k) = (4.13) 
with power spectrum 

(7tk(r)7tk'(r)) = h27rf^6\k + k'). (4.14) 

4.2 Some examples 

Now we will give some examples to show that it is possible to construct states that satisfy 
constraints (4.8) and (4.13). Vacuum state \G) = |0) is obviously the most trivial example 
of a state that satisfies both these constraints. It is also easy to check that coherent states 
with condition (3.3) satisfy both constraints (4.8) and (4.13). 

To give another nontrivial example we look for states of the form \G) = \0) + \^), where 
is an excited state (or a combination of excited states). The simplest example of such a 
state is \G) = |0) + / d^k a{k) a|^|0). But it is easy to check that this state does not work. 
a{k) has to be real to satisfy equation (4.8). Then one can show that a{k) has to vanish to 
satisfy equation (4.13). But 

\G) = |0) + J d^k a{k) al\0) + J d^kid^ks (3{ki)l3{k2) al/jO), (4.15) 

where a{k) and f3{k) are real functions, makes it possible to construct states that satisfy both 
constraints (4.8) and (4.13) and hence produce a scale-invariant power spectrum. Equation 
(4.8) is already satisfied and equation (4.13) leads to 

a{k)a{k') + (4iV - l)/3(fc)/3(A;') = (4.16) 

with = {2'iTfJd^k{[3{k))'^. So aU the states given by equation (4.15) with a{k) and /3{k) 
obeying equation (4.16) produce scale invariant power spectrum. Example of one such state 
is 



\G) = \0) + A^l-^^^^ J d^ke-''' dl\0)+A^ J d^kid^ks 6-^^''^+''^^ 4i4jO)' (4-17) 

^It is possible to avoid the first constraint (4.8) by introducing a new field 7?.j^'^"'(r) — 7?.k('!") — {TZk{T)) and 
defining the power spectrum as {7?.k'^™(r)7?.J^f'"(r)). That will obviously modify the second constraint (4.11). 
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where, A and 7 are real constants and '^^^3'^^ < 1. Renormalizability of the energy- 
momentum tensor of the fluctuations can impose more constraints on these states. Here 
we should note that it is possible to construct infinite number of such states (with even more 
complicated combination of excited states) that obey constraints (4.8) and (4.13) and hence 
can produce a scale-invariant power-spectrum. 

5 Constraints from renormalizability of the energy-momentum tensor 

5.1 Energy momentum tensor for perturbations 

In the longitudinal gauge, the metric has the form 



-(1 + 2$)dt2 + flS^t) [(1 _ 2^)6ij] dx'dxK 



(5.1) 



The Einstein's equation is G^y — S-nGT^y = li^y = 0. To lowest order background obeys 



n 



(0) 



0. The first order Einstein's equations II^^J = give the equations of motion for the 



perturbations. At first order we get $ = Other equations of motion for the perturbations 
are 



+ ms^ + V" 



r,2 



( 47rG02 



+ 



^' =47rG 



V 0, 

We have the good old gauge invariant variable TZ defined by. 



+ 



(5.2) 
(5.3) 

(5.4) 



7^ = If + 



H , 



(5.5) 



Following [35], we can write down the energy- momentum tensor for the perturbations t^i/ 
(2) 

—11^^; in terms of ^ and 6<j) 



too 



SvrG 



+ 



(5.6) 



a'6. 



1 



+ 



3 



42^2 ^ 1(^^)2 _ 1 -2^y^^)2 _ 4^(5^^) _ lv"{4>)i5(j)f + 2V'{^)i^6cP) 

ZD 2 



(5.7) 



Our goal is to find out constraints on the initial state from renormalizability of the energy- 
momentum tensor. So we want to consider the contributions of large-fc fluctuations of tf^^. 
At large- A; (A; >> aH), 

n^r^^-^. (5.8) 

ad) Vk 
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From equation (5.2) we have for large-/c (and for quasi-de Sitter inflation) 

iHea , 



k 



where, e is the slow roll parameter. In this approximation, for 5(j) we obtain. 



IZ- 



H 



(5.9) 



(5.10) 



By inspection, it is very clear that all terms in equation (5.6) containing ^ are suppressed 
by powers of Ha/k compared to terms without ^ . Therefore in large-fc limit we have. 



tij 6ij 



(5.11) 
(5.12) 



5.2 Renormalization of energy-momentum tensor 

Using equation (5.10), we can write down t^i, (for large-/c) in terms of gauge invariant variable 

7^ 



(too) 



H 



(5.13) 



H 



(5.14) 



For the vacuum state, we have the following expressions for the unregularized energy- momentum 
tensor 



(0|too|0) 



1 



(0|t,,|0)«-a%-^ 



dkH^ 
1 



1 + VMl+3^)+r4fe3 



dkH^ 



(5.15) 
(5.16) 



where is the second slow-roll parameter. To obtain the renormalized value of (t^i/), one can 
use any regularization method (for example adiabatic regularization). Detailed discussions 
of adiabatic regularization method can be found in [36, 37]. Adiabatic regularization of {t^y) 
can be done by subtracting adiabatic modes up to order four. 



ren — (0|t^;/ |o) - (o| vlo> 



adi ) 



where, (0|t^i^|0)adi is calculated using the adiabatic mode functions of order four 



(5.17) 



adi-A I 



where. 



w^(^) = k 



1 1 

+ 



u 



adi;i i 



kT^ k^T^' 



(5.18) 



(5.19) 
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Therefore we liave the following expression for (0|t^,y|0) 



I adi 



(0|t00|0)adi 
(0|%|0)adi 



1 



dkH^ 



47r2 



dkH^ 



+ UV finite 
1 



^ + -r^/fcfl + 97?) - -rU 
2k 6 ^ ^ 3 



.4,^3 



(5.20) 

+ UV finite. (5.21) 



Hence, {0\tfj.u\0)ren does not have any ultraviolet divergences. For a general initial state \G) 
we can renormalize energy-momentum tensor by 



(G|V|G) 



{G\i^u\G) - (0|V|0) 



/adi ■ 



(5.22) 



Note that in the right hand side, we do not have {G\t^u\G)adi- But we have {0\t ^,^10) adi- 
It can be understood easily in Minkowski space limit. In Minkowski space, {G\t^^\G)adi = 
{G\tfj,^\G). So, if we had {G\t^^\G)adi in equation (5.22), then {G\t^u\G)ren = for any state 
Clearly that can not be right. 

We will impose the constraint on the initial state \G) that {G\tfj_u\G)ren does not have 
any ultraviolet divergences. That means the state \G) does not introduce any new ultra-violet 
divergences to the energy-momentum tensor. Therefore, to make sure that {G\t^u\G)ren has 
desired UV behavior, we will impose the following constraint on the initial state 



{G\i^,\G) = (0|V|0) +UV finite. 



(5.23) 



If we impose constraint (5.23) on the coherent state ak|C) = C(k)\G) (with C*(— k) = C(k)), 
we find that C(k) goes to zero faster than for large k. One can check that state (4.17) 
also satisfies constraint (5.23). 



6 Bispectrum in the squeezed limit 

In principle one can calculate the three-point function in the momentum space using an 
expression similar to equation (3.11). At leading order in slow-roll parameter, just like the 
coherent state case, we can write 

(G|7t£^(r)7t£;^(r)7t£^(r)|G) ={G\TZ^,{T)TZ^,{r)']Z^,{r)\G) 

- ((G|^k,(T)|G)(G|7tk2(r)7tk3(r)|G) + cyclic perm) . 

(6.1) 

In practice, it is not very easy to compute three-point function for any general state \G). 
Creminelli and Zaldarriaga [38] (see also [39] for a good discussion) used a clever argument 
to determine the three-point function in the local limit (i.e. ki,k2 ^» k^). In the local limit, 
bispectrum has the form 

(7^kl(T)7^k2(r)7^k3(r))fc3«fc,^fc, = {27rf 5^ {ki + k2 + ks)^ fj^lPR{k3)PR{ki). (6.2) 
They found out for any single field inflation model 

/lfi = :|(l-n.) + 0(^)'^ (6.3) 
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It is possible to extend their argument even for any general initial state. We want to compute 
(G|7^ki(''")'^k2(''")^k3(i")|G) after ki,k2 modes have crossed the horizon, ks mode will have 
crossed the horizon in the distant past because ki,k2 ^» k^. Therefore, k^ mode behaves 
classically (see [38]) and contributes to the background metric® 

ds^ = -dt" + a2(t)e-2^^Wdx2, (6.4) 
where, TZb is the contributions from modes far outside the horizon 

7^b(x,t)=/ ^n^e^"^-^. (6.5) 

Therefore, in the local limit [k^ <^ ki ^ k2) we have, 

(G|7^k,(r)7^k,(r)7^k3(T)|G) « (GK^kl(T)7^k,(r))7^,3 7^k3(T)|G). (6.6) 

Where, (7^ki(i")'7^k2(''"))7^k3 — (G'|'^ki(T)'^k2(''")|G)7?,^3 is the two-point function in the per- 
turbed background (6.4). It is clear from the metric (6.4) that long wavelength mode (ne- 
glecting the gradients) is equivalent to the rescaling of the coordinates 

X ^ x' = Ax, (6.7) 

where, A = e~^^. Under this rescaling of the coordinates, it is easy to check that T^k 
transforms as 

7^k ^ A'TlkA = A37^k. (6.8) 

Therefore, 

(G|7tk,(r)7tk,(r)|G)7e,3 = A6(G|%,(r)%,(r)|G) 

= A6(27r)3p;e(fei)<53(ki+k2) 

= A6(27r)3p7e(fei)(A-^+"0'5'(ki + ka)^ 

= (27r)3p7e(fci)(A-i+"0'^'(ki+k2) 

= {27TfPn{ki) (e^^(i-"=)) 5''{ki + ks) 

= {27rfPTi{ki) [1 + 7^B(1 - Us) + ...] 6\ki + ka). (6.9) 

Only thing we have assumed here is that in the superhorizon limit, (G|'^ki (''")^k2 (''")! G^) ^ 
(^^(ki + k2). So far we have treated 7^k3 as a classical field. Now we will promote 7^k3 
to an quantum operator. That leads to 

(G|7tki(r)7tk2(r)^k3(r)|G) « {G\{ni,,{T)n^,{T))n,n^Ar)\G) 

= {27:fPn{kr)6\\i^ + \i2){l-ns) / —^e'''-^ {G\Mr)nU^)\G) 

« {27:fPn{ki)Pn{kz)6\ki + k2)(l - n,) / d^ke'^-^S^k + ks) 

J k<^ki ,k2 

« i27TfPniki)Pn{k3)il - n,)5\-ki + ka). (6.10) 

^Now we are working in the comoving gauge. 
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Therefore, in the squeezed hmit we have, 

(6.11) 

If the state |G) produces almost scale invariant power spectrum then it can not produce large 
This result is actually true for any single field inflation model with a general initial 
state for primordial fluctuations in the limit — 0. For finite k^, the correction to (6.11) 
can be large for some particular states (as reported in [24, 40]). Therefore for finite k^, one 
should be more careful and for a given state a detailed calculation should be performed using 
(6.1). 



7 Conclusions 

In this paper, we have explored the possibility of a general initial state for primordial fluc- 
tuations. Constraints on initial state from current measurements of power spectrum and 
bispectrum are relatively weak and for quasi-de Sitter inflation (or slow roll inflation), a 
large number of states are consistent with the observations. Vacuum state is just one such 
example. Coherent states with constraint (3.3) are also interesting examples of states that 
are consistent with current observations. It is impossible to differentiate between these co- 
herent states and the Bunch-Davies vacuum state just from the power-spectrum and the 
bispectrum. Although we have not considered tensor modes here, similar analysis can be 
done for gravitational waves with similar results. 

Most of the results in this paper only apply to canonical, single-fleld inflationary model 
with equation of state parameter w approximately constant during inflationary epoch. But 
the most signiflcant limitation is the ambiguity in the deflnition of vacuum. We have used 
Bunch-Davies prescription of vacuum throughout the paper. The only justiflcation for it, 
is the fact that for the modes deep inside horizon, space looks Minkowski and hence the 
Bunch-Davies vacuum is a reasonable choice. 

It is clear from our discussion (section 4 and 6) that a large class of general initial 
states can produce scale invariant power spectrum and small just like the Bunch-Davies 
vacuum state. But higher N-point functions can be a very useful tool for probing the initial 
state. The PLANCK satellite will signiflcantly increase the precision of CMB observations. 
Observations made by PLANCK may contain valuable information about the initial state of 
primordial fluctuations and that would provide a window for the physics before inflation. 
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A Computation of three-point function with coherent state 

Using time-dependent perturbation theory we have, 
(C7|7tk,(r)7tk,(r)7tk3(r)|C) 



(C| Te'^-o 



-i r H' jT')dT' 



\C) 



T"0 



(A.l) 



Now all the fields are in the interaction picture and Hint is given by equation (3.10). T and 
T are the time and anti-time ordered product respectively, tq is the conformal time at the 
beginning of inflation and we will take the limit tq — )• — oo.^ Throughout the calculation we 
will assume that kj 7^ 0. We will also take the usual limit r — )• 0. The first term in (A.l) can 
be written using the redefined field (3.9) 



kC\1^..{r)1li^{T)nl^(T)\C) = (C7|7t^,k.(r)^U(r)7t^,k3W|C> 



(A.2) 



+ -(36-2r/) 



(2^^(C|^^,k,(r)7t^,k,(T)7t^,p(r)^^,k3-p(^)|C) + cyclic permj 



1 { f d^v (k3 - P)^ 



(27r)3 kl 



{C\ni^^^{T)ni^^{T)nl^^^^ perm 



7^^j^(r) behaves like the free field, and can be written as 

^lk(T) = ^ [dkniir] + dl^nkir) 

where TZkir) = ^"^^ (l + Using the commutation relation 



■^c,k('^)'4' 



V2 



(A.3) 



(A.4) 



the first term in the right hand side of equation (A.2) can be calculated 

(C|^,^,ki(r)7t^,k,(r)^^,k3(T)|C) = -^(2v^)3C(k3)]2i^e[7^fe3(r)])P^(A;l)53(k^ 

+ ^J(27r)37^fc,(r)7^^^(r)(C|7^^,k,W|C>'^'(k3 + k2) +ki o 



. (A.5) 



In the limit r — 0, (C|7^^ j^(r)|C) = because of the constraint (3.3) and Tlk{T) is purely 
imaginary. Therefore, 

(C|7^^,kl(^)^c,k. W^c,k3(^)|C) = 0. (A.6) 
Last two terms can also be computed and in the limit r — t- we get, 

1 



+ le (pR{k2)PR{h 



(3e - 2rj) (Pr(A;2)Pr(A;i) + cyclic perm) 



-|- cyclic perm 



53(ki + k2+k3). (A.7) 



For large ro, all the exponentials with tq will oscillate. When performing our calculations, we will use the 
average value (i.e. zero) for them. 
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Next we will compute, 



H 



1 



X ( ) [{Re[n'p^{T')]Re[n'^^{T')]C{ps)Cip2)n^^ 



+ p-cyclic) + k-cyclic] — 2e- 



(2vr)S 



1.2 u2 

-J3(k) P^(A:i)Pk(A;2)4;/ + cyclic perm 



'i{ki + k2 + ks 

Where, p = pi + p2 + Ps, k = ki + k2 + ks and we have used the following equation 

/ 1 



(A.8) 



Akr' 



^2 / a{T')k 



(A.9) 



Similarly, 

dT'(CK,(r')7t^,(r)7t^,(r)7t^3(r)|C) = 2e / drV(r'^ 



H 



d'^Pid^Pid^PzS^p) 



1 



X I -2 ) [{Re[l^^^{T')]Re[T^^^{T')]C{p^)C{p2)nk,{T)^^^ 

(2vr)3 



+ p-cyclic) + k-cyclic] + 2e- 



'i{ki + k2 + k-i 

Now we have to compute the last term in equation (3.11) 
{C\n^,{T)\C){C\n^,{T)n^,{T)\C) = 2e [ dr'a^r') 



PR{ki)PR{k2)-j:T^ + cyclic perm 



J,3 

^3 



(A.IO) 



T"0 



x(2i^e[7^;3(rO]i^e[7^;^(r')]C(p3)C(p2)/m[7^*,^(r)7^',^(r')]PR(fc2)<5=^(k2 + k3)<53(ki + pi) 

+p-cyclic). (A. 11) 

Putting all the terms together in equation (3.11), we have 



4e klkl 



1 / kl + kl 

2 + ' kl ) + (A:i + A:2 + A;3) /fef 



+ cyclic perm. (A. 12) 
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